Do mixtures of bosonic and fermionic atoms adiabatically heat up in optical lattices? 
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Mixtures of bosonic and fermionic atoms in optical lattices provide a promising arena to study strongly 
correlated systems. In experiments realizing such mixtures in the quantum degenerate regime the temperature 
is a key parameter. In this work, we investigate the intrinsic heating and cooling effects due to an entropy- 
preserving raising of the optical lattice potential. We analyze this process, identify the generic behavior valid 
for a wide range of parameters, and discuss it quantitatively for the recent experiments with **^Rb and '"'K 
atoms. In the absence of a lattice, we treat the bosons in the Hartree-Fock-Bogoliubov-Popov-approximation, 
including the fermions in a self-consistent mean field interaction. In the presence of the full three-dimensional 
lattice, we use a strong coupling expansion. As a result of the presence of the fermions, the temperature of the 
mixture after the lattice ramp-up is always higher than for the pure bosonic case. This sheds light onto a key 
point in the analysis of recent experiments. 
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Interacting bosonic and fermionic systems play a key role in 
several contexts in physics, quite prominently in the BCS the- 
ory of superconductivity. Systems of dilute atomic gases (in 
optical lattices) offer the perspective of simulating such mix- 
tures or purely bosonic or fermionic systems under extraor- 
dinarily controlled conditions [Ql Bose-Fermi mixtures 
in optical lattices exhibit a rich physical behavior, including 
a wealth of novel phases, charge density waves and super- 
solids ijstl. Recent experiments have succeeded in preparing 
such mixtures in optical lattices ['^J, notably the realization of 
a stable bosonic '^^Rb and fermionic ''^°K mixture in three- 
dimensional optical lattices IHS]. 

To achieve realizations of such strongly correlated systems 
in the quantum-degenerate regime, very low temperatures 
have to be reached. This is not only a difficult prescription but 
also, while thermometry methods in the absence of a lattice 
are established, it is not entirely clear how to measure the tem- 
perature in its presence. Indeed, following recent experiments 
with cold bosonic atoms, an intriguing and fruitful contro- 
versy |70 has arisen concerning the general question relevant 
to experiments with ultracold atoms in optical lattices; How 
cold, after all, is the system in the optical lattice expected to 
be? For Bose-Fermi mixtures, this question is even harder to 
answer as the additional degrees of freedom leave more room 
for different explanations. Interactions between bosons and 
fermions result in an effectively reduced repulsion between 
bosons, independent of the sign of the Bose-Fermi interaction. 
Hence, one might well expect an increase in coherence as 
compared to the purely bosonic case. Quite surprisingly, how- 
ever, the opposite effect (as measured in terms of the visibil- 
ity of the quasi-momentum distribution) was observed iHI^. 
The theoretical work Ref. Ql, based on numerical quantum 
Monte Carlo and density-matrix-renormalization-group sim- 
ulations of one-dimensional systems, points towards the pos- 



sibility that this might actually be due to a finite temperature 
effect. 

In this work, we discuss the thermodynamics of adiabatic 
loading of harmonically trapped Bose-Fermi mixtures into op- 
tical lattices. During the adiabatic loading procedure, the en- 
tropy remains constant and leads to intrinsic cooling or heat- 
ing processes. We argue that one should expect a significant 
adiabatic heating of the mixture, not to be confused with ex- 
perimental imperfections such as parametric heating. This is 
by no means a marginal effect, as we will quantitatively clar- 
ify. This resulting temperature determines the physics of the 
strongly correlated system once the optical lattice is present. 

We identify the generic behavior and discuss it on the basis 
of the values corresponding to the experiment described in 
Ref. 101 ■ More precisely, the presence of fermions leads 
either to a more distinct heating of the mixture or a less 
distinct cooling. We study in detail the behavior of these 
adiabatic heating and cooling effects in the inhomogeneous 
system — complementing results for purely (non-)interacting 
bosonic (Jgtl) ifiol [Till and non-interacting fermionic lll2ll 
systems — and analyse and flesh out the specific role of the 
fermions in this adiabatic process. 

Trapped Bose-Fermi mixture without optical lattice. - 
Subsequently, we will discuss the thermodynamics of the 
Bose-Fermi mixture in an isotropic harmonic trap in the 
absence of an optical lattice. We will insist on being close 
to an experimental situation in our description, and take the 
full three-dimensional situation into account. We start from 
the grand-canonical Hamiltonian 
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where we denoted the bosonic (fermionic) field operators by 
$ (^P), the interaction amplitudes gsB, gpB are related to 
the respective scattering lengths as gsB = '^T^f^CLBB/mB, 
gpB = 2Tih?aFB{'m'B + mi?)/ (mi? ms), the free part of the 
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FIG. 1: Densities qb = uo^^/Nb (red), gp = maj^^/Np (green) 
and entropy densities (Jb/f = Sb/f<4ioI^^* (blue/cyan, scaled by 
10'' for clarity), atot = era + ctf (black), in units of the bosonic 
harmonic oscillator length aho- The shown data sets correspond to 
the experiment in Ref. Q] and a temperature of 95 nK. The mixture 
consists of 10^ *^Rb atoms and no (solid), Nf = Q.QZNb (dashed), 
and Nf = Q-QINb (dotted) ''"K atoms. 



bosonic Hamiltonian is given by 



IniB 



and accordingly for hp. We thus restrict ourselves to isotropic 
traps, taking geometrical averages of the trapping frequencies 
in the actual experiment 

For the bosonic sector, we invoke the standard local density 
Hartree-Fock-Bogoliubov-Popov (HFBP) approximation, 
which is a self-consistent mean field scheme that has proven 
applicable to a wide temperature regime, see Refs. ifm [itIi 
and references therein. The inter-species interaction is treated 
in the self-consistent mean-field approximation (see, e.g., 
Refs. dlMl) 

where we defined the fermionic density m and the total 
bosonic density ?i = tiq + ?it composed of the condensate 
and non-condensate density hq, ht, respectively. This 
yields the following set of coupled equations: (i) The finite 
temperature Gross-Pitaevskii equation in the Thomas-Fermi 
approximation (which may be safely applied for the high 
number of atoms considered), governing the condensate 
density 



uq = max < 



^J.B - Vb - QFBm 
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where the chemical potential is fixed by the given total number 
of bosons, Nb = Nq + Nt — JdruQ + Jdrnr- (ii) The 



thermal density of bosons (fc^T =!//?) 



riT 



dp 
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where the Bogoliubov amplitudes are given by 



2e 

and the quasi-particle spectrum reads 
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2m_B 



Vb - f^B + 2gBBn + gFBm ] ~ Qbb^o 



Finally, (iii) the fermionic density in local density 
approximation 



(27r)-'Mp 
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Vf - fJ-F + gFBfl, 



where the chemical potential is fixed by the given total number 
of fermions Nf = /dr to. 

For given temperature T and particle numbers Nb, Nf, we 
solve (i)-(iii) self-consistently in the following way: Starting 
with no interaction between bosons and fermions, gpB = 
and riT = 0, we (a) compute no and /is by solving (i) under 
the particle number restriction, (b) obtain nr from (ii), (c) 
iterate (a) and (b) until convergence, (d) solve (iii), which 
yields m and /ip, (f) iterate (a)-(d) until convergence. 

After convergence, we are equipped with the energies e, 5, 
and can compute the entropy of the mixture, 

S/kB = J [sb{p, r) + SFip, r)] , 
where the individual contributions read 



^B(p,r).-i^-logfl-e-''^(-) 



SFip,r) = 



PSip,r) 



log 1 + e 



This is the expression forming the starting point of the 
analysis in the absence of the lattice. In Fig. [T] we show 
the obtained results for the parameters of the experiments in 
Ref. ^ for different ratios Nf/Nb- The critical temperature 
for Rose condensation is w 205 nK for all three cases 
Nf/Nb ~ 0,0.03,0.07. In this experiment no thermal 
cloud was discernible, corresponding to a EEC fraction of 
at least 80% and a initial temperature below 95 nK. The 
bosonic entropy is highest at the condensate boundary, where 
the density of the thermal cloud has its maximum. In turn. 
Si? is highest in the center of the trap. We can see that 
the bosonic contribution, sb, to the total entropy remains 
basically unaltered by the presence of the fermions, their main 
contribution to S stemming from s f itself. 

Trapped mixture in deep optical lattices. - To describe the 
system in the presence of the lattice, we use the single-band 
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Bose-Fermi-Hubbard Hamiltonian lUl |2ll] H = J + hi, 
where 

hi = ^ni{hi - 1) + Vh.riii - iifh, - fif Wj. 

Here, the operator bi (fi) annihilates a boson (fermion) at 
site i and fii = Sj^S,;, = ///i. J accounts for the 
tunneling of atoms from one site to neighboring sites, U, V are 
the intra-, respectively inter-species on-site interactions, and 
/i - ' = /ig jp — ViWQ on-site chemical potentials controlling 
the particle number via /i^ jp and accounting for the harmonic 
confinement V^, which is approximately the same for both 
species. For deep lattices the tunneling may be treated as a 
perturbation in Jb/f- Up to second order and within local 
density approximation (assuming that the trapping potentials 
are the same at neighboring sites), the free energy is found 
to be = - \og{Z)/(3 = -E.log(2.0//3 -iY..{Jlh + 
Jlh)lzl where 6, = T.Zn'=,<^' + ^K'""' ^ 
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and the unperturbed on-site energies and corresponding 
partition functions are given by e"'™ = Vnin — 1) — [ifn — 

/if m + Vnm, Zi = Er=oNp(-/3er'°) + exp(-/3e"'^)]. 
The structure of a sum over on-site terms originates from 
the local-density approximation and helps to speed up the 
computation substantially: All sites with the same value Vi 
yield the same contribution to the free energy. They may thus 
be grouped together and the contribution be calculated only 
once per group as opposed to for every single site. Starting 
from the above expression for the free energy, we calculate 
the chemical potentials for given particle numbers Ng /f by 
numerically solving Ng/p = —dF/dfiB/F^ where the right 
hand side is obtained by numerically differentiating the free 
energy with respect to the chemical potentials. Similarly, by 
differentiating with respect to /3, we then compute the entropy 
of the mixture in the lattice: S/ks = P^dF/dp. 

Discussion. - We are now in the position to assess the 
situation when ramping up the lattice in an adiabatic process. 
Fig. |2] shows the entropy as a function of temperature in a 
system of 10'^ ^^Rb-atoms in a three-dimensional trap. In this 
figure, we use the experimental parameters of Ref. but 
the findings are valid for a wide range of parameters. Both 
the entropy without and in the presence of the optical lattice 
is depicted for the purely bosonic case and a small admixture 
of fermions iVj^/TYs = 0.03,0.07. We see that generally, 
for fixed lattice depth, below a certain temperature, the 
adiabatic ramp-up gives rise to an adiabatic cooling, whereas 
above this threshold temperature we find an adiabatic heating. 
Both above and below this temperature, the presence of the 
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FIG. 2: Entropy as a function of temperature for the parameters of 
Ref. 0]. The mixture consists of 10^ **^Rb atoms and A^f = 
(solid), O.OSA^B (dashed), Np = 0.07Nb (dotted) '"'K atoms in a 
lattice of various depths (blue: ISEr, red: 30Er, black: no lat- 
tice). The inset shows the same at a larger scale, including the en- 
tropy for the purely bosonic case obtained from analytical expres- 
sions (green) valid below the critical temperature (dashed) and at 
ultra-low T (solid), see footnote jT4ll . 



fermions results in a higher final temperature as compared to 
the same situations with bosons only. This is most dramatic 
at initial temperatures for which without fermions the bosons 
are adiabatically cooled and in the mixture adiabatic heating 
occurs, e.g., for an initial temperature of 90 nK and a final 
lattice depth of l^En the temperature is w 40 nK higher in 
the presence of Np = Q.07Nb fermions, corresponding to an 
increase of « 67 %, see Fig. [3] This affects the contrast of the 
interference pattern analyzed in those experiments. This 
behavior is generic, valid in particular for both experiments 
of Refs. mia] as well as for experiments performed in an 
isotropic and shallower trap ifisll . where the initial temperature 
was always below the threshold: For any initial temperature, 
the entropy without the lattice is always much higher in the 
presence of fermions, even for the relatively small admixture 
of '"^K atoms as in Fig.|2] While below the threshold adiabatic 
cooling occurs, this effect is lessened compared to the purely 
bosonic case. Above the threshold and in the presence of the 
lattice, the entropy including fermions is higher, thus reducing 
the heating effect. This can however not compensate for the 
high initial difference of entropies, see Fig. [3] 

Note that the influence of fermions is most distinguished in 
the absence of the lattice. This is plausible when considering 
the form of the unperturbed free energy in the presence of the 
lattice: e"'° = Un{n — 1) — fifn and e"'^ = Un{n — 1) — 
/if 71 + Vn are different only by an alteration of a definition 
of the bosonic chemicals potentials (the total number of 
bosons is the same with and without fermions), leading for 
low temperatures to approximately the same expression for 
the entropy. Taking a closer look at the situation including 
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FIG. 3: Difference AT/i™ = 2;f„^ - Tyf„ between the final tem- 
perature in the lattice with and without fermions as a function of the 
initial temperature without lattice. Parameters are as in Fig. |2l For 
any initial temperature, the presence of the fermions leads to higher 
final temperatures as compared to the purely bosonic case. Solid 
lines depict the final temperature TJif„ without fermions (right scale). 



the lattice, we see that, at low temperatures, more fermions 
lead to a lower entropy - the inter-species attraction reducing 
the mobility of the atoms and thus reducing the number 
of possible micro-states. In turn, at higher temperatures 
interactions become less important and the entropy increases 
with the number of fermions. While this can not compensate 
the initial difference in entropies, it however reduces the 
heating effect for higher initial temperatures, see Fig. [3] 
This effect may be observed in the currently available 
experiments: At a fixed lattice depth, the difference between 
the situation with and without fermions should first increase, 
reach a maximum and finally decrease with increasing initial 



temperature. 

Summary and outlook. - In this work, we have 
quantitatively explored the adiabatic cooling and heating 
effects that are to be expected in experiments with Bose-Fermi 
mixtures in optical lattices, crucial when reaching a strongly 
correlated system. On intuitive grounds, one could have 
suspected that the features observed in experiments were 
entirely due to a shift of the bosonic Mott lobes in the 
presence of fermions, the presence of fermions effectively 
altering the local chemical potential. This is indeed the 
case, but predicts an increase of coherence i^, the opposite 
of which was observed in experiments. We have seen that 
under the parameters as used in present experiments, the 
resulting temperature is much larger than expected from 
thermometry based on measurements before the ramp-up of 
the lattice. Methods to assess the temperature of samples 
within deep optical lattices would clearly be a breakthrough 
for any studies on quantum gases in lattices. Promising 
ideas are, e.g., the detailed characterization of the shell 
structure of local densities ifTsll . A link to the expected 
visibility from our analysis is provided by Ref. Isd. This 
analysis applies to a one-dimensional situation, yet for the 
visiblility, it is expected to give a clear guideline: It is 
seen how the bosonic visibility decreases as the temperature 
increases. A clear-cut quantitative analytical analysis of 
the quasi-momentum distribution at finite temperature is still 
lacking and poses — even for purely bosonic systems — an 
exciting challenge and constitutes a test-bed for theories 
developed in the condensed matter context. It is the hope 
that the present work can significantly contribute to the 
clarification of the intriguing discussion on the interpretation 
of observed data and on the available theoretical models for 
ultracold mixtures of bosonic and fermionic atoms in optical 
lattices. 
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